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We compute the force acting on a free, static electric charge outside a uniform dielectric sphere. 
We view this force as a self-interaction force, and compute it by applying the Lorentz force di- 
rectly to the charge's electric field at its location. We regularize the divergent bare force using 
two different methods: A direct application of the Quinn-Wald Comparison Axiom, and Mode-Sum 
Regularization. 



PACS numbers: 41.20.-q, 41.60.-m 



I. INTRODUCTION 



Electric charges which are placed in an inhomogeneous 
ponderable medium undergo self interaction. The sim- 
plest case is that of a static electric charge in an inho- 
mogeneous dielectric. The self interaction of the charge 
results in a self force (in other contexts this self force is 
also known as a radiation reaction force), which acts to 
accelerate the charge. In the static problem, then, one 
can ask the following question: What is the external, 
possibly non-electric, force which needs to be exerted on 
the charge to keep it static when a non-uniform dielectric 
is present? 

In this paper we study this question for a very sim- 
ple case. Specifically, we find the self force on a point- 
like electric charge e outside a uniform dielectric sphere. 
(By an inhomogeneous dielectric we here mean the dis- 
continuity of the dielectric constant at the surface of the 
sphere.) The origin of the force on the charge in this case 
is simple: The charge polarizes the dielectric at order e. 
The induced electric field then back-reacts on the origi- 
nal charge, and this interaction then is at order e 2 . In 
fact, one can compute the force on the charge following 
this simple physical picture. However, one can employ 
a different picture, in which the force on the charge is 
construed as a self force. The charge interacts with its 
own field, and the latter is distorted by the presence of 
the dielectric sphere. In this picture the force is com- 
puted locally, using only the fields at the location of the 
charge. The local approach has many merits. Specifi- 
cally, the computation of the near field is much simpler, 
and one does not have to compute additional quantities 
such as the far field or the sphere's polarization. (The 
fields at great distances of course contribute to the force 
on the charge, but for the local approach only through 
boundary conditions.) A similar approach was used in 
Ref. Q to compute the radiated power of synchrotron 
radiation using only the near field. The difficulty in the 
local approach arises from the well known fact that the 
field of a point charge diverges when the evaluation point 
for the field coincides with the field's source. In fact, this 
happens already for a static charge in empty (and flat) 
spacetime. Let the charge e be located on the z-axis in 
spherical coordinates. Decomposing the scalar potential 



$ into Legendre polynomials, one finds that 



00 p 
— rt 



$(r,tf) = e^Tfr^(costf), 



where r < (r> ) is the smaller (greater) of the r values of 
the source's location ro and the evaluation point. The 
(self) force acting on the charge is then given by the 
(average of the two one-sided) gradients of the potential. 
Specifically, 
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which clearly diverges. (The derivation of the last equa- 
tion is given below.) In this illustration of the problem, of 
course, it is clear that the regularized, physical self force 
vanishes: The force obviously cannot depend on where 
we choose to put the origin of our coordinate system. 
(Also, we have ample observational evidence that static 
isolated charges remain static.) 

There is a long history of works on the self force. (For 
reviews see, e.g., ||.) Recently, the analogous problem 
of calculation of self forces in curved spacetimes (also for 
the gravitational case where the self interaction pushes 
a body with finite mass off a geodesic) has gained much 
interest ]3|, |j . In this paper we shall make use of some of 
the techniques, which have been developed for self inter- 
action in curved spacetime, for the problem of interest. 
(Interestingly, there is a close link between electromag- 
netism in static gravitational fields and electromagnetism 
in matter. As is well known ||, Maxwell's equations 
in vacuum in static curved spacetime can be written as 
Maxwell's equations in flat spacetime with an effective 
non-uniform dielectric.) Specifically, we shall make use 
of the Quinn-Wald Comparison Axiom || and Mode-Sum 
Regularization jj] in order to extract the physical, finite 
piece of the self force. 

The organization of this paper is as follows. In Section 
II we solve for the scalar potential, and obtain the modes 
of the bare force. This is, in fact, a standard exercise in 
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FIG. 1: An electric charge e at a distance ro from the center 
of a sphere of radius R. The sphere has a dielectric constant 
e = 1 + eo, and outside the sphere the dielectric constant is 
unity. Region I is for r < R, region II for R < r < ro, and 
region III for r > tq. 



electromagnetism J7[. Then, in Section HI we regular- 
ize the self force using two different approaches, and in 
Section IV we discuss the properties of our result. 



II. DERIVATION OF THE BARE FORCE 

Consider a static electric charge e in vacuum at radius 
ro, outside an insulated sphere of radius R of uniform 
dielectric constant e = 1 + eo, where e > 0. Notice, 
that eo is not the permittivity of free space, but rather 
e = 47%, where %e is the electric susceptibility. We 
place the charge e on the z-axis without loss of generality. 
This configuration is plotted in Figure [jj 

Maxwell's equation in matter is 



V D = Airp 



(1) 



where D = eE is the displacement field, E is the electric 
field, and p is the density of free charges. 

We assume that the distribution of the dielectric is 
spherically symmetric (although non-uniform). Specifi- 
cally, we take e = e(r). (Despite the uniformity of the 
sphere, the dielectric constant throughout space depends 
on r: it suffers a step-function discontinuity at the sur- 
face of the sphere.) In the usual spherical coordinates 
Eq. (nj) becomes 
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where $ is the scalar potential. We next decompose 
Eq. (|J) into Legendre polynomials. That is, $(r, d) = 
£ € ^(r)P*(oos0) and p(r,§) = ^±^1^(21 + 
l)P((cos{l). The radial equation then becomes 
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The boundary conditions for this equation are that <j/ is 
continuous everywhere (which includes regularity at the 
origin and at infinity), but d r (j> 1 is discontinuous at r = R 
and at r = ro. Specifically, these latter two conditions 
are that 

lim e{R-a)d r cj>\R-a)= lim e(R + a)d r <j) l {R + a) 

a— *0+ a— >0+ 

(which comes from the continuity of the normal compo- 
nent of the displacement field at the surface of disconti- 
nuity), and 

lim [d r <j>*{r + n) - $y (r„ - a)] = -{21 + 1)4, 

(7 — >0+ Tq 

[which comes from integration of Eq. (^) across r = ro , 
and using the continuity of 4> e and e there (the only dis- 
continuity of e is at r — R)]. 

The radial functions <f> 1 then satisfy 

{Air r < R (region I) 

B e r e + Cir- 1 - 1 R<r < r (region II) ,(4) 
Dtr~ l ~ x r > ro (region II) 

where the coefficients Ai,Bi,Ce and Di are found from 
the boundary conditions. We find that 
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such that the scalar potential <!> is given by 
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Here, 



$ vac - E w+T+lTo ^Pt^Pz (cos 0) r>R 



.(9) 
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is the potential in the absence of a dielectric sphere. 

The bare force f bare is found by f barc = -eV<I>, evalu- 
ated at the location of the charge at r = tq and § = 0. 
From symmetry, it is clear that any force is radial. We 
compute, then, the radial component of the force only. 
Differentiating Eq. (||) and using Eq. (flo|), we find that 
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^vac to construct the fields F^. Applying the Compari- 
son Axiom, we find that the self force is given by 



fr = -E 
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where /* = __| lim CT ^ 0+ [<9 r </(r + er) + d T 4> l (r Q - a)]. 
Clearly, Eq. ([ll]) diverges. This comes as no surprise, as 
we have already mentioned that this divergence occures 
already for a charge in empty space. In the next section 
we shall extract the physical, finite part of this infinite 
bare force. 



III. REGULARIZATION OF THE BARE FORCE 

In order to regularize the bare force (|TT|) , we make di- 
rect use of the Quinn-Wald comparison axiom, for which 
plausible arguments were given. The Comparison Axiom 
states the following (see M for more details): Consider 
two points, P and P , each lying on time-like world lines 
in possibly different spacetimes which contain Maxwell 
fields F^y and F^ sourced by particles of charge e on the 
world lines. If the four- accelerations of the world lines 
at P and P have the same magnitude, and if we iden- 
tify the neighborhoods of P and P via the exponential 
map such that the four-velocities and four- accelerations 
are identified via Riemann normal coordinates, then the 
difference between the electromagnetic forces f^ and f^ 
is given by the limit x — > of the Lorentz force associ- 
ated with the difference of the two fields averaged over a 
sphere at geodesic distance x from the world line at P , 
i.e., 

- U = lim e (f»u - F^) u v . (12) 

x— >0 \ I x 

Here, we identify the "tilde" spacetime as that of a 
globally empty spacetime. Obviously, / M = 0. We em- 
phasize that this axiom assumes a nearly trivial form for 
the case of interest: The local neighborhood of the par- 
ticle in question and of a similar particle in a (globally- 
)empty spacetime are identical. (It is only the far-away 
properties of spacetime - as represented by different di- 
electric constants - which are different for the two space- 
times.) Another remark is that we do not need to average 
here over directions, as the forces in our case are direc- 
tion independent. Consider now Eq. for the potential. 
Outside the dielectric sphere the potential $ contains the 
vacuum potential 'I'vac and a correction A<i>. We next use 



2-Fi being the hypergeometric function. We were unable 
to find this result in the literature. (In view of the vast- 
ness of the literature on classical electromagnetism, our 
search in the literature is naturally incomplete.) 

Before we analyze the properties of this result, let us 
derive it using a second method. Specifically, we use 
Mode-Sum Regularization. (Note, that Mode-Sum Reg- 
ularization is based on the Quinn-Wald result for the self 
force in curved spacetime, the latter being a consequence 
of the Comparison Axiom. In that sense, these two 
methods are not entirely independent. Here, however, 
we make direct use of the Comparison Axiom, which is 
necessary but not sufficient in order to derive the Quinn- 
Wald result.) Mode-Sum Regularization is described in 
Refs. [|. In Mode-S um regularization one finds two reg- 
ularization functions, and d M . The regularized self 
force is given by 



/M = E(4 bale -^)-^ ( 15 ) 

1=0 

where is a finite valued function and h e ^ has the gen- 
eral form h E r — a r (£ + i) + b r + c r (£ + i) -1 . One only 
needs the local properties of spacetime in order to deter- 
mine these functions. As locally the charge is in empty 
space (it is removed from the dielectric sphere), it is clear 
that the regularization functions hL and d M would be the 
same as in a globally-empty spacetime. Indeed, it is easy 
to find the limit as t — > oo of the modes of the bare force. 
The modes of the radial component of the bare force 
( |i"l| ) approach e 2 /(2rg) as I — > oo. As hL must have the 
same asymptotic structure (as i — > oo) as this implies 
that h l T = e 2 /(2r l ), identically the same as in (globally- 
)empty spacetime, in agreement with the previous rea- 
soning. We similarly expect the function d r to vanish, 
as it does in a globally-empty flat spacetime. We justify 
this expectation a posteriori by demonstrating that this 
leads to the same expression as we received by using the 
Comparison Axiom. It then follows that the regularized 

self force is given by f r = J2Ur haTC ~ e2 /(2 r o)]i which 

I 

agrees with Eq. (|l3"|). 
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For any finite value of eo the force is smaller in magni- 
tude than in the case of a conducting sphere (|l6|) . This 
behavior is shown in Fig. |]A, which plots the self force 
as a function of eo for fixed tq. It can be seen that as 
eo - * oo, the full expression approaches the saturation 
value of the conducting sphere. 

At very large distances (ro S> R), the self force be- 
comes 
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FIG. 2: The self force on a free charge outside a dielectric 
sphere. The charge e is taken to be that of an electron, and 
the radius of the sphere is R — 1cm. Upper panel (A): The 
self force as a function of eo, for ro = 2cm. Lower panel (B): 
The self force as a function of ro (in cm), for eo = 10.7. 



IV. PROPERTIES OF THE RESULT 

We found that the self force on the charge e is given 
by Eq. This is an attractive force, as indeed is 

expected. (The charge e polarizes the sphere such that 
there is an excess of oppositely-charged induced charge 
on the sphere closer to the free charge. Hence the polar- 
ization charge acts to attract the free charge.) We can 
check our result in the limiting case of infinite dielectric, 
eo — * oo, which corresponds to the case of an uncharged, 
insulated, conducting sphere. In that limit our result 
becomes 
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which is indeed the known result for an uncharged, insu- 
lated, conducting sphere §. The opposite extreme case 
is the limit as e — * 0. Linearizing our result in e , we 
find that 



f r = - ^ 
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which vanishes linearly with eo as eo —* 0. 



which drops off like r§ . This behavior can be seen from 
Fig. ||B, which displays the self force as a function of ro 
for fixed eo. We can check the validity of this limit by 
deriving Eq. (^J|) using the alternative picture [|| . When 
ro ^S> R, we can treat the field of the charge e to the 
leading order in l/r*o as a constant over the sphere. Let 
us take for simplicity the charge e = sgn(e)|e| to be on the 
positive z-axis. The electric field due to the free charge is 
Eo = — sgn(e)|e|rQ~ 2 z, and the polarization of the sphere 
then is just a constant inside the sphere, and is given by 
P = [3/(47]-)] [e /(3 + e )]E . The dipole moment p can 
be obtained by a volume integral over P. One finds then 
that p = [eo/(3 + e )]i? 3 E . The electric field E at tqZ is 
found by E = [3z(p • i) — p]/Vg, or E = — 2prg~ 3 sgn(e)z, 
where p = |p| = eoi? 3 rg" 2 |e|/(3 + e ). The force on the 
charge e is simply f = eE = — 2jpr n ~ 3 |e|z, which is equal 
to the leading order term of Eq. ( J18| ) . 

When r approaches R the self force grows rapidly, 
and in the limit diverges. This is indeed expected: in this 
limit one has a point charge near a semi-infinite dielectric. 
The solution for the force is a classic image problem [0, 
which obviously diverges in the coincidence limit of the 
charge and its image. This divergence happens already 
in the case of the conducting sphere, as is evident from 
Eq. (16). In fact, we find that the self force diverges 
whenever the free charge is locally at a region with non- 
zero gradient of the dielectric constant. 

We note that the magnitude of this self force is not 
extremely small for realistic parameters. Take the charge 
e to be that of an electron of mass m e , and the dielectric 
sphere to be made of Silicon, for which eo = 10.7 at room 
temperature and pressure, and take the sphere to be of 
radius 1cm. In the gravitational field of the Earth, with 
gravitational acceleration of 980cm/sec 2 , the self force 
equals the weight of the electron when r = 13.2cm 
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